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Abstract 

The low temperature solution of the exact master equation for an oscillator coupled 
to a linear passive heat bath is known to give rise to serious divergences. We now 
show that, even in the high temperature regime, problems also exist, notably the 
fact that the density matrix is not necessarily positive. 
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1 Introduction 



In a previous publication [1], we presented a general solution of the exact 
master equation for an oscillator coupled to a linear passive heat bath. This 
was achieved by solving the generalized quantum Langevin equation for the 
initial value problem [2] which enabled us to obtain explicit expressions for the 
time-dependent coefficients occuring in the exact master equation. Mext, we 
illustrated the general solution with the construction of explicit expressions 
for the probability distributions for two examples: wave packet spreading and 
"Schrodinger cat" state for a free particle. This enabled us to discover that the 
low temperature solution gave rise to divergences. The purpose of the present 
paper is to show that problems also exist in the high temperature regime. 
Thus, we are motivated to extend our previous analysis to obtain explicit 
expressions for the Wigner function and the elements of the density matrix, 
and with the aid of these results to critically examine the equation and its 
solution. 

At the outset we must point out that the exact master equation has drawbacks 
that seriously limit its utility. For the most part these are related to the 
fact that the initial state is necessarily one in which the particle and heat 
bath are uncoupled, while the subsequent motion is described by the fully 
coupled system. The most serious of these is that the equation (and therefore 
its solution) has an irreparable divergence, irreparable in the sense that it 
persists in a cut-off model with finite bath relaxation time [1] . We believe this 
divergence is related to a well known phenomenon of quantum field theory: 
the Hilbert space spanned by the eigenf unctions of the coupled Hamiltonian 
is orthogonal to that of the uncoupled Hamiltonian [3] . It can be shown that 



2 



this is the case for the system of coupled oscillators that is the microscopic 
basis for the derivation of the exact master equation and its solution [4,5]. 

The exception is the high temperature limit, where the divergence does not 
appear since, by convention, one there neglects the zero-point oscillations of 
the bath. We therefore confine our discussion to that limit. However, even in 
that case there are further difficulties. The first of these is that the initial state 
breaks the translational invariance of the system. This is manifest in an initial 
squeeze centered at the origin that, for a free particle, results in a drift of the 
mean position toward the origin. For this reason in the examples we consider 
only the case of a free particle, where this difficulty is most apparent. 

A second difficulty with the equation in the high temperature limit is that 
for short times the diagonal elements of the density matrix, which have the 
physical interpretation of probabilities, are not necessarily positive. This is 
a well known difficulty with the weak coupling master equation, solved by a 
further time average to bring it to the so-called Lindblad form that guarantees 
positivity [6]. It was supposed that the exact master equation, with its time- 
dependent coefficients, would be immune but, as we shall show explicitly, in 
the high temperature limit this difficulty remains. The situation, therefore, 
is that there are concerns that seriously limit the utility of the exact master 
equation. The irreparable divergence can only be avoided with the neglect of 
zero-point oscillations in the high temperature limit, and even there problems 
remain with the failure of translational invariance and positivity. 

Despite these problems, we argue that the solution of the exact master equa- 
tion can be useful. In particular, we show that when the initial particle temper- 
ature is taken to be equal to that of the bath, we obtain results that for short 
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times are identical with those obtained from a calculation which considers en- 
tanglement at all times [7]. Finally, an important result is the demonstration 
that measures of decoherence based on the Wigner function are identical with 
those based on the off-diagonal elements of the density matrix. We argue that 
the Wigner function, being everywhere real, is the preferable description. 



2 General solution in high temperature limit 



The exact master equation, in terms of the Wigner function W(q,p,t) for the 
complete system of oscillator plus heat bath, may be expressed in the form 

dW 1 dW „ 2/ s dW 

-xt = P^— + ™i (tjQ^r- 

ot m oq op 

dvW d 2 W d 2 W 

+ 2T(t)^- + hmT(t)h(t)^- + m)f(t)^- p , (1) 

and we note the presence of four time-dependent coefficients for which we 
obtained explicit expressions [1]. 

The general solution of the exact master equation is most simply expressed 
in terms of the Wigner characteristic function (the Fourier transform of the 
Wigner function). 

W(Q,P;t) = JdqJ dpe- i{qP+pQ)/h W{q,p;t). (2) 

Expressed in terms of this function, the general solution given in Eq. (4.15) 
of [1] takes the simple form, 



(X 2 ) P 2 + m (XX + XX) QP + m 2 (X 2 ) Q 2 
W(Q, P- 1) = exp{--^ \ 

xW(mGQ + GP,m 2 GQ + mGP;0). (3) 
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In this expression X(t) is the fluctuating position operator, defined in Eq. 
(2.16) of [1], while G(t) is the Green function, defined in general to be 

G(t 2 - h) = hx{h), x{t 2 )]8(t 2 - h), (4) 
in 

in which x{t) is the time-dependent Heisenberg coordinate operator and 6 is 
the Heaviside function. 

It was shown in [1] (Sec. V.A.2) that the quadratic expectations (X 2 ), (XX + XX 
and (x 2 ^ are logarithmically divergent, which has the effect that the solution 
(3) vanishes identically for non-zero positive times. There, too, it was shown 
that the divergence is irreparable, in the sense that it persists in a model with 
a high frequency cut-off. As was stated in the Introduction, it appears that 
this is a reflection of the well established fact that the states of the uncou- 
pled system (the initial state) are orthogonal to those of the coupled system 
(states for later times) [4,5]. This fact, surprising from the point of view of 
the quantum mechanics of finite systems, is a feature only of systems with an 
infinite number of degrees of freedom. In our case it is the heat bath that must 
be infinite. In any event, the upshot is that, while the derivation of the exact 
master equation is formally correct, there is no useful result. The exception 
is the high temperature limit, where by convention one neglects zero-point 
oscillations. We therefore restrict our further discussion to this limit. Indeed, 
it is this limit (or better said: approximation) that has been considered in all 
previous discussions of the exact master equation. 

In considering the result in the high temperature limit, we restrict our discus- 
sion to the Ohmic model of coupling to the heat bath, which is adequate for 
our purposes. This model corresponds to Newtonian friction, with retarding 
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force proportional to the velocity. In addition we consider only the case of a 
free particle, moving in the absence of any external potential. Choosing the 
friction constant to be 777.7, the Green function then takes the simple form: 



1 - e~ 7 * 

G(t) = — . (5) 

7777 



In the high temperature limit the moments appearing in the solution can then 
be written in the form: 



1 

X 2 ) = 2 mi kT J dt'G 2 (t') 




kT 

[ 27 t_(i- e -*)(3- e -*)], 



777.7 2 



XX + XX) = 2m^kTG 2 {t) 

2kT _ 
■(1-e 



7 t\2 



7777 

t 



X 2 ) = 2m 1 kT J dt'G 2 {t 
kT 







— (1-e-n (6) 
777 



We should emphasize that in these expressions T is the temperature of the 
bath. In the initial state the bath and the particle are uncoupled, so in the 
general solution (3) W(Q, P; 0), the initial Wigner characteristic function, may 
be chosen to have any form consistent with a single particle not coupled to 
a bath. In the examples we consider two possibilities: the particle initially at 
zero temperature and at a temperature equal to that of the bath. 
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3 Example: Gaussian wave packet 



3. 1 Initial particle temperature zero 

We consider first an initial state corresponding to a Gaussian wave packet, 
with wave function of the form 

^Hp^e X p{-<^24 (7) 

This is a stationary wave packet, centered at x with mean square width a 2 . It 
is a minimum uncertainty wave packet, since the mean square momentum has 
the minimum value H 2 /4a 2 . The corresponding Wigner characteristic function 
is 

Q nw*, , Q 



W{Q,P;0)= f alxe- iPx/h (f){x-^-,0)(f)*{x+^-,0) 
J 2 2 



Q 2 a 2 P 2 x P 
= 6Xp{ -8^-^-^ } - (8) 



This initial state is what is termed a pure state. The condition for a pure state 
is usually stated in terms of the density matrix: Tr{p 2 } = 1. In terms of the 
Wigner characteristic function, this condition becomes 

oo oo 

— JoIqJ dP\w(Q,P)\ 2 = l. (9) 

— oo — oo 

It is easy to show that a state is a pure state if and only if it can be associated 
to a wave function. Such a pure state necessarily corresponds to a particle at 
zero temperature; a finite-temperature state is a mixed state. 
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Putting the expression (8) for the initial Wigner characteristic function in the 
general solution (3) we see that 



.xo(m 2 GQ + mGP) , , . 

W(Q, P; t) = exp{ — % }, (10) 



where we have introduced 



A^ = (x 2 ) + a 2 m 2 G 2 + ^- 

kT h 2 (l - e^) 2 

ft r [2 7 t - (1 - e -7*)(3 - e"*)] + a 2 e~ 2 ^ + 1 ' 



m / y 2 ' 4m 2 7 2 a 2 ' 



,11 + 11 . .. h 2 mGG 

2 4cH 

= K(l _ e -7*)2 _ m7C7 2 e -2 7 t + ft 2 (e^-e 2 ^) 

7 4m 7 cr i 



A^=m 2 (x 2 ) + a 2 m'G 2 + ^^- 

= mkT(l - e~ 2 ^) + m 2 7 Ve- 2 ^ + (11) 



Here we have used the superscript "(0)" to indicate that the particle is initially 
at zero temperature before suddenly being coupled to the bath at temperature 
T. 

We now form the Wigner function with the inverse Fourier transform 

W(q,p; t) = JdQj dPe^ p+ ^ h W(Q, P; t). (12) 

With the expression (10) for W(Q, P; t), the integration is a standard Gaussian 
integral [1], with the result 

W(q,p;t) = W (0) (q - mGx ,p - m 2 Gx ;t), (13) 
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where 



1 j(%2 _ 94(0) , j(0) 2 

W {0) (q, P ;t) = . 1 exp{- 229 Jf% W+ (0 i llP }(14) 

2^a?M^ - 4f 2(4° 1 ) 4? - 4f) 



is the Wigner function corresponding to a single minimum uncertainty wave 
packet initially at the origin. 

The Wigner distribution (13) corresponds to a Gaussian distribution in phase 
space, centered at the point ({x(t)) , (p(t))), where 

(x(t)) = mGx = x e~ Jt , 

(p(t)) = m 2 Gx = -m7x e" 7t . (15) 

The width of the distribution is characterized by the coefficients , which 
have the interpretation 



Aff = A^ = ^(f))-<*(f)r, 
A^ = Ap 2 = ( P 2 (t))-( P (t)} 2 , 



4? = \ (x(t)p(t) + p(t)x(t)) - (x(t)) (p(t)) . (16) 

The breaking of translational invariance is evident, since the center of the 
distribution drifts to the origin, which is a special point. We get more insight 
into this phenomenon by forming the limit as t goes to zero through positive 
values, 

uil n+\ 1 J (g-^o) 2 2a 2 (p + m 7 g) 2 

W( g ,p;0 ) = ~ h eM ^ jr 2 }• (17) 

Here we see that the particle has received an impulse —m^q proportional to 
the displacement and directed toward the origin. This is exactly the action that 
produces a squeezed state [8,9]. Indeed, as one can readily verify, the Wigner 
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function (17) corresponds to a pure state with associated wave function 

^•^(^^{-^-^ < 18 > 

Thus, the state immediately after the coupling to the bath is still a zero 
temperature pure state. However, as a result of the squeeze toward the origin, it 
is no longer a minimum uncertainty state. It still corresponds to a wave packet 
centered at x with mean square width a 2 , but the mean square momentum 
uncertainty is now Ap 2 = h 2 /4cr 2 + m 2/ ~f 2 a 2 , greater than the minimum value. 
The result of the squeeze is that the distribution (13) drifts toward the origin. 
Less obvious perhaps, the same phenomenon appears in a shrinking width of 
the probability distribution at short times. At longer times the width expands 
due to thermal spreading 

3.2 Initial particle temperature equal to that of the bath 

As we have noted, the minimum uncertainty wave packet (7) considered above 
corresponds to a particle at temperature zero. A perhaps more realistic initial 
state would be that for a particle initially at temperature T, equal to that of 
the bath. In [1] we showed that the initial state at temperature T is obtained 
from that at temperature zero by the replacement 

W(Q, P- 0) -> exp{-^^}iy(Q, P- 0). (19) 

This finite temperature state is a mixed state, no longer satisfying the condi- 
tion (9) for a pure state. It cannot be associated to a wave function. Making 
this replacement with the zero temperature form (8) and repeating the steps 
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leading to the result (13), we obtain 

W{q,p;t) = W {T) (q-mGx ,p-m 2 Gx ;t), 



(20) 



where now the Wigner function corresponding to a single Gaussian wave 
packet initially at the origin takes the form 

1 A^q 2 - 2A^ T \p + A^p 2 

w (9 ' p;t) = 2^AS-A^ eM - M 2«^-V) K(21> 

in which 



A^=A ( S+mkTG 2 

2kT . 1 _ 2 ^ 2 (l-e-^) 2 
= — (jt - 1 + e"* + (J 2 e- 2 ^ + ^ 2 2 / , 

4? =4? +m 2 kTGG 

kT h 2 e~^{\ — e~^\ 

= (1 - e -7t) _ m7 a 2 e- 27t + / 2 J , 

fc2 -2<yt 

= m/cT + m 2 7 Ve" 2 ^ + (22) 

Note that the motion of the center of the distribution is the same as that 
given in Eq. (15) for the case of initial particle temperature zero. The only 
change is an increase in the width of the Gaussian distribution due to thermal 
spreading, the widths being given by the expressions (16) with A^} in place 
of Ajjj. We again get more insight into this by forming the limit as t goes to 
zero through positive values, 



W(q, p; 0+) = exp{ — h2 — }. (23) 

nhJl + 4& 2a 2(^ + mkT) 



Here we see that that there is the same initial squeeze as in Eq. (17 ), but the 
squeeze acts on a state with mean square momentum uncertainty increased 
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by the mean square thermal momentum, mkT. The result is that immediately 
after the squeeze the mean square momentum uncertainty is Ap 2 = h 2 /4cr 2 + 
mkT+m 2 ^ 2 a 2 . However, the initial Ax 2 = a 2 and is unaffected. by the squeeze. 

4 Example: Pair of Gaussian wave packets 

4-1 Initial particle temperature zero 

The wave function for an initial "Schrodinger cat" state consisting of two 
separated minimum uncertainty Gaussian wave packets has the form 

«*■<» = (8 ^)i/4 (1 + e -,V^) 1 / 2 (^{- i ^>+^P{- if 4^ ! >)'( 24 ) 

The initial Wigner characteristic function is then given by 

Therefore, from the general solution (3) we can write, 

x (cos (gggg^gg)j + cosh (^+^)^ 6) 

Here the A$ are again given by (11). 

The inverse Fourier transform (12) again involves only standard Gaussian 
integrals. We thus easily see that the Wigner function for a wave-packet pair 
is 
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W( q , P ; t) = 2{1 + e . dV8 ^ \ q --P-—) 

+2e- AW{t) W i0 \q,p) cos$ (0) (g,P : t)}, (27) 

where is the Wigner function (14) for a single minimum uncertainty wave 
packet at the origin and we have introduced 



A<°>(f) 



f a(0) h 2 G 2 \j A (°) h 2 m 2 G 2 \ f a(0) h 2 mGG\2 



4<H 



4 (0) 4 (0) .(0)2 g 2 

^11 ^22 ^12 



$ w>p.*) = 4 (o), ( o) A m 4^2- ^ 

■^11 ^22 ~~ ^12 

The first two terms in brackets in (27) correspond to a pair of Gaussian wave 
packets of the form (13), centered initially at x = ±d/2 and propagating 
independently. We can call these the direct terms. The third term is an in- 
terference term. The direct terms are a maximum at their centers, which are 
drifting toward the origin but for d ^> a will for short times be well away from 
the origin. The interference term is a maximum at the origin and initially its 
peak height is exactly twice that of either of the direct terms. (However, as 
one can readily verify, the area under the interference term is independent of 
time and a factor e ~ d2 / Sa2 smaller than the area under the direct terms.) The 
factor e~ A(0) ^\ equal to the ratio of the peak height of the interference term to 
twice the peak height of the either of the direct terms, is taken as a measure 
of the interference in phase space. In Fig. 1 we show as a function of 
and there we see that for <C 1, is linear in t. Using the expressions (11) 
for the quantities Afi and (5) for G we can readily show that for short times 

A< >(f)^7«, (29) 
A th 
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where A t h is the mean de Broglie wavelength, 

Ath = T i=. (30) 



For a separation large compared with the thermal de Broglie wave length, this 
corresponds to a decay of the peak height of the interference term on a time 
scale short compared with 

This is the frequently appearing "decoherence time" [7]. 

It is of interest to form the probability distribution, which in general is given 
by 

00 

P(x;t)= J dpW(q,p;t) 

—00 

00 

= 7r^ [ dPe ixP/h W(0,P;t). (32) 

—00 

Using the expression (26) for the Wigner characteristic function, this can be 
written 



r m 2 G 2 d 2 , / nW , HGdx . . . 

+2a(t) exp{— _^-}pW(x) cos 7^(0) )' ( 33 ) 

where 

^ ( ° ) (^) = ^i^exp{--^ y } (34) 



is the probability distribution for a single minimum uncertainty wave packet 
initially centered a the origin and we have introduced the attenuation coeffi- 
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cient, 



(X 2 ) d 2 
8a 2 A { u 



a(t) = expj-^^}. (35) 



The attenuation coefficient is a measure of the interference as it would be 
observed in the probability distribution. We should emphasize that the proba- 
bility distribution can be directly measured. This is to be contrasted with the 
Wigner distribution in phase space, which is not directly observable. 

The integrated probability in the interference term is equal to e~ d2 / 8 °" 2 , inde- 
pendent of time. For d ^> a this is negligibly small. But it is the amplitude 
of the interference fringes that is observed and measured by the attenuation 
coefficient (35). 

For short times, <C 1, using the expressions (6) for (X 2 ) and (11) for A$ , 
we see that the attenuation coefficient takes the form 



a ^- eXp{ - 3r d [^ + (WA) 2 ] } - (36) 



This initially falls very slowly from unity, but after a time 2ma 2 /H (the time 
for the mean square width Ax 2 to double) will decay with a characteristic time 
3rd, where is the characteristic time (31) for the decay of the interference 
term in the Wigner function. 

With no coupling to the bath the interference pattern in the probability dis- 
tribution would not appear until the direct terms begin to overlap due to wave 
packet spreading. Taking 7 — > in the expression (11) for the mean square 
width, we see that this would be a time of order t mix = 2mad/h <C r d . What 
this means is that neither the rapid disappearance of the interference term 
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in the Wigner distribution, nor the corresponding rapid decay of the attenu- 
ation coefficient, can could be directly observed. All that can be seen is the 
non-appearance of the interference term in the probability distribution. 



4-2 Initial particle temperature equal to that of the bath 



The corresponding results for the particle initially at the temperature T of 
the bath are obtained using the prescription (19). The result is 



mGd m 2 Gd, 



„, m , mGd m 2 Gd. 



+2e" A(T) {t) W {T) (q, p) cos ¥ T \q,p:t)}, (37) 



where 



( a(T) _ h 2 G 2 \( A T ) _ h 2 m 2 G 2 \ _ ( a(T) _ h 2 mGG \2 j2 
A(T)( f \ _ I 71 !! 4a 2 ^22 4 a 2 > 1^12 4a 2 > " 

{ ) A (T) A {T) - A {T)2 8a 2 

A U A 22 A 12 

{T) _ (GA%> - mGA$)q + (mGA™ - GA$)p fid 

1 j ~ 4? 4? - A { S 2 i^" 1 ' 

This expression for the Wigner function has the same form as that found above 
for the case where the initial particle temperature is zero. The only difference 
is that 4ifc — ¥ ^ e area unc ^ er the interference term is still independent 
of time and smaller than the area under the direct terms by the same factor 
e -d 2 /8cr 2 mam diff erence i s that the peak height of the interference term 
is much reduced. The ratio of the peak height of the interference term to twice 
that of the direct terms is given by the factor e~ A(T) ^ but, as we see in Fig. 
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1, no longer vanishes at t — 0. Instead, we now find for <C 1, 

A< r )(t) 9* x / 2 - . (39) 



For separation large compared with both the thermal de Broglie wavelength 
and the slit width, this will be a large number and the peak height of the 
interference term will be correspondingly small. Thus, we see that there is 
no time scale for the disappearance of the interference term, it is already 
negligibly small at t — 0. We might say that this is just the point. What 
has occurred is that an incoherent superposition of pure states has wiped 
out the interference term. In the present case this is the result of the initial 
preparation of the state. In the previous case of an initial pure state, this 
occurs as a result of the "warming up" of the particle in a time r^. In either 
case the interference term is gone long before its effect could be observed in 
the probability distribution. 

The probability distribution is of the same form (33), but with A u (t) replac- 
ing A^i(t). The attenuation coefficient is now given by 



((X 2 ) + mk TG 2 )d 2 
8cr 2 4i 



a(t) = exp{- " ; 7 ' } (40) 



For yt 1 this becomes 

■M^-I^i < 41 > 



These results for short times are identical with those obtained from an exact 
calculation assuming entanglement of the particle with the bath at all times 
[7]). Since the coupling to the bath, as measured by the decay rate 7, does 
not appear they can also be obtained from a calculation without dissipation 
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using only the familiar formulas of elementary quantum mechanics [10]. Thus 
we see that, while the exact master equation gives here correct results, it does 
so only for times so short that coupling to the bath can be neglected. 

5 Density matrix 

The density matrix is a Hermitian operator, p, defined such that for any wave 
function ip corresponding to a state of the particle in the Hilbert space (ip, pip) 
is the probability that the system is in the state. From the notion of probability, 
we see immediately the the density matrix must be a positive definite operator. 
Consider the eigenf unctions <p a and the corresponding eigenvalues p a of the 
density matrix, 

P<Pa = Pa<Pa- (42) 

Clearly, the eigenvalue p a is the probability that the system is in state <p a and 
is therefore necessarily positive. We assume the density matrix is normalized, 
so that 

Efa = l- (43) 
a 

This is a sum of positive terms equal to unity. It follows immediately that 
Tr{p 2 } = E^<l- (44) 

a 

Moreover, the equality holds if and only if there is exactly one eigenstate 
with probability one, all others having probability zero. This is exactly the 
condition (9) for a pure state. 
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In coordinate space the elements of the density matrix are given in terms of 
the Wigner characteristic function by the formula: 



oo ; 

(x\p\x')= J d V e i{x - x ' )p/h W(^^ lP ) 

— oo 

oo 

2^E I dPel(X+X ' )P/2h W(x'-x,P). 



The normalization condition (43) becomes 



Tr{p} = J dx(x\p\x) = #(0,0) = 1. 



(45) 



(46) 



That is, the Wigner characteristic function corresponding to a normalized 
density matrix must have the value unity at the origin. 

Next consider 



oo oo 



Tr{p 2 }= J dx J dx'(x\p\x') (x'\p\x 

— oo — oo 

00 oo 



(47) 



-oo — oo 



where we have used that fact that W(—Q, —P) = W(Q,P)*. Thus, the con- 
dition (44) takes the form 



oo oo 



! r dQ f dP W(Q,P) 2 < 1, 



2nh 



(48) 



-oo — oo 



with the equality holding if and only if the density matrix is that of a pure 
state. This is the condition (9) for a pure state. 
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5. 1 Gaussian wave packet 



Consider now the case of an initial state corresponding to a single minimum 
uncertainty wave packet, for which the Wigner characteristic function at time 
t is given in Eq. (10). The evaluation of the formula (45) for the matrix element 
of the density operator involves a single Gaussian integral. The result can be 
written 

m 2 G 

(x \p(t) \ x') = exp{i — - — x (x — x')}R(x — mGx , x' — mGxo; t), (49) 

where we have introduced 

p/ 1 , M<^-AT)(x-x') 2 + (x + x') 2 , .. A^-x^) 
R(x,x ;t) = , exp{-- ^ h in [ 

which is the matrix element corresponding to a single wave packet initially 
centered at the origin. Using the expressions (5) for the Green function and 
(11) for Ajk(t), we see that viewed in the xx' plane the absolute value of the 
density matrix, \(x \p(t) \ x')\, consists of a single Gaussian peak. This peak is 
initially centered at x = x' = x and circularly symmetric with width a. In 
the course of time the center of the peak drifts to the origin while stretching 
in the direction x — x'. 

For short times, the density matrix is not necessarily positive. The simplest 
way to see this is to form Tr{p 2 }, using the above expression for the density 
matrix element or evaluating directly the expression (47). In either case we 
find 

Tr{p 2 } = j (51) 

z V /1 ll /1 22 A 12 
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In Fig. 2 we show a plot over a short initial time interval of this result, 
evaluated using the expressions (11) for the with parameters A t h = 4cr, 
kT = hftrf. There we see clearly that the condition (44) is violated for short 
times. It follows that for such times there must be negative eigenvalues; the 
density matrix is not positive-definite. We gain some insight into this phe- 
nomenon if we use the expressions (11) to expand 

r 2 



TrV 2 } = l + (l-4^ r ) 7 t + ---. (52) 



Hh 

From this we see that Tr{p 2 } will increase above unity whenever A th > 2a. 
The failure of positivity occurs for a sharply peaked wave packet. 

To see explicitly this failure of positivity, it will be sufficient to take xq = 0. We 
then form the expectation of the density matrix with a wave function ip(x) oc 
d<t>(x, + )/dx, where <p(x, + ) is the state immediately after the coupling to 
the bath (given in Eq. (18) with xo — 0). Explicitly, 

^=^w eM - ( ^ +i ^ )x2Y (53) 

With this, 



oo oo 

(V,p(*)^)= / dx j dx'^*(x) (x\p(t)\x')^(x') 



—oo — oo 



2\l - ±(A {0) A {0) - A {0)2 )] 



[(1 + An)(l + ^(A^^-AT ) ) + !MVt4?)! ]3/2 

This expectation will be negative whenever A(A^A^ — A^ 2 )/^ 2 < 1, that is 
exactly for those times when Tr{p 2 } given by (51) is greater than one. 

We have seen that general solution of the exact master equation results in a 
density matrix that can have negative eigenvalues. How can it be that an exact 



(54) 
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result can lead to such an unphysical consequence? The answer, of course, is 
that we been able to give a non-trivial meaning to the solution only in the 
high temperature limit, which involves a serious approximation (neglect of 
zero-point oscillations). It is this approximate result that fails the positivity 
test. 



5.2 Wave packet pair 



For the wave packet pair the Wigner characteristic function at time t is given 
by (26). With this the formula (45) for the matrix element of the density 
operator becomes a standard Gaussian integral. The result is 



(x\p(t)\x>) 1 



2(l + e- rf2 /8^) 
x [exp{iL^-(x — x')}R(x — mG^-, x — mG^-) 

Z Z Z 

d , mGd . mGd 

exp{-tL-(x-x)}R(x + — ^— ,x + — ^— ) 

+e- A(0> exp{tM^(x + x')}R(x - K^,x' + K^) 
z z z 

+e~ A(0> exp{-iM^(x + x')}R(x + K^, x' - K^)}, (55) 
_ z z 

where R(x,x';t) is is the matrix element (50) corresponding to a single wave 
packet initially centered at the origin and to shorten the expressions we have 
introduced 

h*(mGA<$-GA<8) m*G h{mGA^-GA^) 



Viewed in the xx' plane, the density matrix element (55) will shows four peaks 
exactly of the form (49) of single wave packets; two diagonal peaks centered 
at x — x' — ±mGd/2 and two off-diagonal peaks centered at x — —x' = 
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±Kd/2. The off-diagonal peaks are multiplied by the same factor e~ A(0) ^ 
that multiplies the interference term in the Wigner function (37). The result 
is that under the condition d ^> A t h, while the peaks slowly broaden and drift 
toward the origin, the off diagonal peaks rapidly shrink in amplitude. Thus, the 
two descriptions of the disappearance of the interference term, that in terms 
of the Wigner function and that in terms of the density matrix element, are 
entirely equivalent. In this connection note that the diagonal density matrix 
element is the probability distribution (33). The description in terms of the 
Wigner function has the advantage of being real, so somewhat simpler. Neither 
the Wigner distribution in phase space nor the density matrix element in xx' 
space is directly observable, so either description is theoretical. 

6 Conclusion 

We have exhibited a number of significant failures of the exact master equa- 
tion: Due to the irreparable divergence the solution strictly does not exist 
except in the high temperature limit. Even in that approximation the solution 
breaks translational invariance. Finally, the solution corresponds to a den- 
sity matrix that can have negative expected values. But we must remember 
that the derivation of the equation as well as its solution are formally correct. 
Rather, the failures can all be traced back to the assumption of an unen- 
tangeled initial state. The picture of a particle suddenly clamped to a bath 
corresponds to no physically realizable operation. Even when the initial par- 
ticle temperature is adjusted to that of the bath it receives a violent squeeze 
toward the origin. From this follows the "unphysical" consequences: correctly 
described by the exact master equation. Unfortunately, the final assessment 
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is that the exact master equation is of very limited utility. 

In fact, in an exact description the particle must be coupled to the bath (en- 
tangled) at all times. Think of the example of the blackbody radiation field, 
which is coupled to the particle through its charge that can in no way be 
switched on or off. In physical applications the coupling to the radiation field 
is sufficiently weak that a description in terms of the familiar weak coupling 
master equation (of Lindblad form) is appropriate, but the effects of renor- 
malization and zero-point oscillations are nevertheless present at all times. 
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Fig. 1. The functions A^(t) and A^it) multiplied by \ 2 t Jd 2 plotted vs. jt. The 
parameters chosen are A t h = 4cr and kT = Shj. 
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Fig. 2. Tr{p 2 } for an initial minimum uncertainty wave packet. The parameters 
chosen are Ath = 4c, kT = 5hj, the same as for Fig. 1. Note that at short times 
Tr{p 2 } > 1. 
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